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1  Overview 

Reactions  involving  N  and  O  atoms  dominate  the  energetics  of  the  reactive  air  flow  around 
spacecraft  when  reentering  the  atmosphere  in  the  hypersonic  flight  regime.  One  of  the  im¬ 
portant  determinants  for  the  chemistry  in  such  flows  are  the  thermal  rate  coefficients  for 
reactive  processes.  In  the  present  project,  rate  coefficients  are  determined  for  reactions  in¬ 
volving  0(^P)  and  NO(^n).  For  this,  a  potential  energy  surface  (PES)  for  the  ground  state 
of  the  NO2  molecule  is  constructed  based  on  high-level  ab  initio  calculations  and  interpolated 
using  the  reproducible  kernel  Hilbert  space  (RKHS)  method  and  Legendre  polynomials.  The 
global  PES  of  NO2  in  the  ground  state  is  constructed  by  smoothly  connecting  the  surfaces 
of  the  grids  of  various  channels  around  the  equilibrium  NO2  geometry  by  a  distance-based 
switching  function.  The  rate  coefficients  are  calculated  from  Monte  Carlo  sampling.  The 
results  indicate  that  at  high  temperatures  primarily  the  lowest  electronic  state  is  relevant  and 
zero-point  effects  are  not  relevant.  The  rate  coefficient  for  O2— production  (the  “0102-l-N” 
channel)  at  20000  K  becomes  comparable  (to  within  a  factor  of  around  three)  to  the  rate 
coefficient  of  the  oxygen  exchange  channel  for  the  same  temperature.  The  computational 
approach  outlined  and  pursued  is  generic  for  small-molecule  reactions  and  can  be  applied  to 
other  reactions  relevant  to  the  hypersonic  flight  regime. 

2  Introduction 

The  computation  of  effective  rate  coefficients  at  arbitrary  thermodynamic  conditions  is  an 
essential  aspect  in  modeling,  understanding  and  characterizing  the  chemistry  in  the  hypersonic 
flight  regime.  Direct  experimentation  at  the  relevant  conditions  (temperatures  of  up  to  20000 
K;  nonequilibrium  flow)  is  usually  difficult  or  even  impossible.  Hence,  an  in  silico  approach 
holds  much  promise  in  providing  the  necessary  input  data  for  modeling  the  chemical  networks. 

Nitric  dioxide  (NO2)  is  one  of  the  molecules  which  plays  an  important  role  as  a  model 
system  in  atmospheric  chemistry,  as  a  smog  constituent  and  in  combustion  processes.  [1,  2,  3] 
While  most  of  these  processes  occur  close  to  thermal  equilibrium  and  at  moderate  temper¬ 
atures,  reactions  involving  NO2  have  also  attracted  interest  under  more  extreme  conditions. 
This  is  the  case  for  the  hypersonic  flight  regime  of  spacecraft  reentering  the  atmosphere. 
The  chemistry  at  the  surface  of  such  vehicles  typically  involves  highly  non-equilibrium  con¬ 
ditions  with  vibrational  and  rotational  temperatures  reaching  several  thousand  Kelvin.  [4] 
The  gas-phase  and  surface  reactions  and  energy  transfer  at  these  temperatures  are  essen¬ 
tially  uncharacterized  and  the  experimental  methodologies  capable  of  probing  them  are  not 
well  established.  Under  such  circumstances,  validated  computational  investigations  become 
a  valuable  complementary  tool. 

In  the  present  work,  the  formation  of  NO2  from  NO(^n)  and  0(^P)  and  its  subsequent 
dissociation  into  oxygen-exchange  and  O2  formation  channels  is  investigated.  Characterizing 
the  kinetics  and  dynamics  of  a  chemical  reaction  from  computation  usually  requires  a  potential 
energy  surface  for  the  interactions  involved.  Within  the  framework  of  Molecular  Dynamics 
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(MD)  simulations,  this  can  either  be  an  ab  initio/semiempirical  MD  (AIMD)  approach,  a 
model  potential  energy  surface  (PES)  or  a  parametrized,  multidimensional  PES  fitted  to 
electronic  structure  calculations,  experimental  observables,  or  a  combination  of  the  two.  [5] 
Using  AIMD  is  usually  only  possible  for  small  systems  and  sufficiently  short  simulation  times 
due  to  the  considerable  computational  cost. [6,  7,  8]  A  more  efficient  approach  is  to  use  separate 
grids  for  the  different  regions  of  interest  and  join  them  smoothly.  In  the  asymptotic  regions, 
that  is,  reactants  or  products,  the  dynamics  is  carried  out  on  a  single  and  simpler  PES 
leading  to  a  faster  simulation.  In  the  intermediate  regions,  the  global  surface  is  represented 
as  a  weighted  sum  of  the  contributions  of  all  individual  PESs.  This  approach  has  been 
successfully  applied  to  connect  different  force  fields.  [9]  We  will  pursue  a  similar  strategy  in 
the  present  work  to  build  a  global  PES  for  NO2. 

The  aim  of  the  present  work  is  the  characterization  of  the  reaction  between  0(^P)  and 
NO(^n)  at  higher  temperatures  relevant  to  the  hypersonic  flight  regime  of  reentering  space¬ 
crafts.  At  a  more  fundamental  level,  we  are  particularly  interested  in  establishing  an  efficient 
computational  framework  capable  of  probing  small  molecular  systems  for  a  wide  range  of  pa¬ 
rameters.  The  general  approach  is;  (1)  the  construction  of  PESs  for  the  ground  state  of  the 
NO2  molecule  based  on  high-level  ab  initio  calculations  and  their  fitting  with  a  reproducing 
kernel  Hilbert  space  (RKHS)  method  combined  with  Legendre  polynomials;  (2)  quasi-classical 
trajectory  calculations  to  study  the  adiabatic  reaction  dynamics,  and  (3)  calculation  of  the 
rate  coefficients  for  the  different  exit  channels  using  Monte  Carlo  method.  The  results  of 
the  simulations  will  be  compared  with  previous  computational[10,  11,  12]  and  experimental 
works[13]  in  order  to  validate  the  computational  approach.  In  the  context  of  the  hypersonic 
flight  regime,  the  O2  production  channel  (NO-l-0— 7-02-l-N)  is  of  particular  interest  since  the 
thermal  dissociation  threshold  for  O2  is  lower  than  that  for  N2  in  air.  A  central  question  is 
how  translational  excitation  of  the  reactants  and  vibrational  excitations  of  the  NO  affect  the 
rate  coefficient  of  O2  formation. 

3  Methodology 

3.1  Ah  initio  calculations 

The  present  calculations  use  Dunning’s  standard  correlation-consistent  polarized  quadruple 
(cc-pVQZ)  basis  set. [14,  15]  For  a  meaningful  description  of  the  geometries  and  energetics 
of  the  NO2  molecule,  its  dissociation  energetics  and  the  conical  intersections  between  the 
ground  and  the  low-lying  excited  states,  a  multiconfigurational  method  is  required.  For  this, 
the  complete  active  space  self-consistent  field  (CASSCF)  method  is  employed  which  takes  into 
account  static  correlation.  For  recovering  the  dynamic  correlation,  e.g.  during  bond  break¬ 
ing  and  formation,  multireference  configuration  interaction  (MRCI)  is  used.  The  Davidson 
correction  [16,  17]  is  included  to  estimate  the  effects  of  higher  excitations  (MRCI-I-Q).  All 
calculations  are  performed  with  the  MOLPRO  program  ver.  2012.1. [18] 

Figure  1  gives  a  schematic  representation  of  the  general  relevant  stages  of  the  reaction 
of  interest.  Initially,  the  0(^P)  atom  approaches  the  NO(^n)  molecule.  Depending  on  the 
impact  parameter  b,  the  quantum  state  of  the  NO  (vibrational  and  rotational)  and  the  collision 
energy,  the  NO2  molecule  is  formed.  In  the  absence  of  additional  collisional  partners  (zero 
pressure)  the  NO2  will  decay  as  it  was  formed  with  an  energy  above  its  dissociation  limit.  The 
following  possibilities  exist  for  the  products;  NO2  can  (1)  dissociate  into  the  same  incoming 
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Figure  1;  Relevant  stages  in  the  NOl+02  reaction.  The  energies  were  derived  from  the  ab 
initio  calculations  of  this  work. 

constituents  (NOl+02— 7-N01+02),  (2)  dissociate  into  the  same  constituents  but  with  the 
oxygen  atoms  exchanged  (NOl+02— )-N02+01),  (3)  form  dioxygen  (O2+N)  or  (4)  atomize 
into  the  individual  atoms  (N+01+02).  The  latter  is  not  further  considered  in  this  work. 

In  the  presence  of  additional  collisional  partners,  the  NO2  can  loose  its  excess  energy  and 
can  stabilize.  In  the  high-pressure  limit  this  takes  place  for  all  trajectories  which  enter  the 
potential  well  of  NO2  and  thereby  it  is  straightforward  to  determine  the  high-pressure  limit 
rate  coefficient  for  NO2  formation. 

In  order  to  describe  the  different  channels,  three  grids  in  Jacobi  coordinates  are  employed 
to  systematically  cover  the  relevant  regions  of  the  PES  in  the  NO2  system  (Figure  1).  The 
first  grid  (for  the  NOI-l-02  asymptote)  includes  the  coordinates  ri  (the  N-01  distance),  Ri 
(the  distance  from  the  center  of  mass  of  NOl  to  the  oxygen  atom  02),  and  ai  (the  angle 
formed  by  01,  the  center  of  mass  of  NOl  and  02  (Figure  1)).  With  this  grid,  the  entrance 
(NOl-l-02)  and  N02-formation  channels  can  be  described.  For  the  oxygen  exchange  channel 
(N02-01),  a  permutation  of  the  oxygens  in  the  first  grid  is  performed.  As  a  result,  the  second 
grid  is  defined  by  r2,  R2  and  02.  For  the  third  grid  (0102-l-N),  is  the  0-0  distance,  R3  is 
the  distance  from  the  center  of  mass  of  the  O2  molecule  to  the  nitrogen  atom  and  03,  is  the 
angle  between  02,  the  center  of  mass  of  the  O2  molecule  and  the  nitrogen  atom  (Figure  1). 

For  the  ri  and  r3  coordinates,  12  points  {nr  =  12)  from  1.7  to  2.55  ao  were  selected.  As 
for  the  Ri  and  R3  coordinates,  29  points  {ur  =  29)  between  1.85  and  18.9  ao  were  selected. 
The  and  03  angles  (in  degrees)  are  from  an  11-point  (n^  =  11)  Gauss-Legendre  quadrature 
(a  =  11.98,  27.49,  43.10,  58.73,  74.36,  90.0,  105.64,  121.28,  136.9,  152.5,  168.02).  In  total, 
3828  points  were  calculated  for  each  grid. 
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3.2  Representation  of  the  PESs 

In  order  to  obtain  the  PESs  of  the  NO2  molecule  from  the  ab  initio  energies,  a  reproducing 
kernel  Hilbert  space  (RKHS)  ansatz[19]  for  the  distances  (r  and  R)  is  combined  with  Legendre 
polynomials  for  the  angles  (a)  (see  Figure  1  for  the  Jacobi  coordinates).  First,  for  each 
combination  {R,a)  in  the  grid,  a  set  of  functions  VR^aii")  is  obtained  by  using  a  RKHS.  In  a 
second  step,  since  the  ab  initio  calculations  were  carried  out  at  Gauss-Legendre  points,  it  is 
straightforward  to  expand  the  2d-PES  W{R,  a),  in  terms  of  Legendre  polynomials  Px{cosa) 
as  follows, 

10 

W{R,a)  =  ^  fx{R)Px{cosa),  (I) 

A=0 

where  fxiR)  are  radial  coefficients  and  A  is  their  order.  Again,  an  RKHS  interpolation  is  used 
along  the  R  coordinate.  For  the  global  PES  for  the  system,  the  local  PESs  determined  for 
the  three  grids  corresponding  to  the  relevant  channels  of  the  reaction  (Figure  1)  are  joined. 
The  global  PES  is  thus  the  weighted  sum  of  the  three  asymptotic  PESs: 

3 

=  '^Wj{x)Vj{x)  (2) 

i=i 

where  x  =  (ri2,  ri3,  r23)  are  the  three  interatomic  distances  of  the  NO2  molecule  (Figure 
I),  and  Ki(x),  V2{x.)  and  V3(x)  are  the  potential  energies  for  the  three  regions  for  which  the 
distances  ri2  (N-01),  ri3  (N-02)  and  r23  (01-02)  are  small,  and  the  remaining  distances  vary 
from  small  to  large,  respectively.  The  weights  rcj(x)  are  coordinate-dependent  and  result  from 
the  normalization  of  the  following  weight  functions. 


Wjo{xj)  =  exp(— Xj/Ai?), 

(3) 

(4) 

E  Wioixi) 

l=i 


where  AR  =  0.05  uq.  A  similar  weighting  function  was  successfully  used  in  Multisurface 
Adiabatic  Reactive  Molecular  Dynamics  (MS-ARMD)  which  uses  the  energy  as  the  control 
parameter  in  connecting  multiple  force  fields.  [9] 

3.3  Molecular  Dynamics  Simulations 

Armed  with  a  globally  valid  PES,  the  dynamical  evolution  of  the  system  is  followed  by  propa¬ 
gating  Hamilton’s  equations  of  motions  subject  to  initial  conditions.  The  equations  of  motion 
are  integrated  numerically  using  the  Velocity- Verlet  algorithm.  [20]  Since  the  total  linear  and 
angular  momenta  are  conserved  during  the  dynamics,  it  is  convenient  to  select  the  center  of 
mass  of  the  NO2  molecule  as  the  origin  of  the  coordinate  system.  In  order  to  capture  the 
fastest  nuclear  movements  in  the  system,  a  time  step  of  2  a.u.  (a.u.  means  atomic  unit  of 
time)  (4.84  x  10“^  fs)  was  used  for  the  propagation.  Trajectories  were  started  from  an  initial 
separation  of  (min(15  ao,  impact  parameter  -|-  5  ao))  and  followed  until  the  fragments  had 
separated  to  the  same  distances.  If  this  criterion  was  not  met,  integration  continued  for  a 
maximum  time  of  the  estimated  interactionless  time  of  flight  through  the  “reaction-zone”  plus 
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an  additional  100  ps,  during  which  most  of  the  NO2  molecules  decayed  except  for  very  few 
low-energy  (low  and  low  rovibrational  state  of  initial  NO)  trajectories. 

Suitable  initial  conditions  for  NO  were  generated  from  a  WKB[21]-quantized  periodic  or¬ 
bit  of  the  corresponding  rotating  Morse  oscillator  for  given  vibrational  v  and  rotational  j 
quantum  numbers.  [22]  The  symmetry  axis  of  the  NO  molecule  and  the  axis  of  its  rotation, 
which  are  orthogonal  to  each  other,  was  also  randomly  orientated  as  was  the  NO  angular  mo¬ 
mentum.  Sampling  of  the  initial  rovibrational  quantum  states  of  NO,  the  impact  parameter 
(6)  and  the  collision  energy  {Ec)  is  discussed  in  the  next  section. 

For  characterizing  the  reaction  at  a  given  temperature,  4  different  cases  (I  to  IV)  are  consid¬ 
ered: 

(I)  corresponds  to  NO2— formation  in  the  high-pressure  limit.  The  NO2  molecule  stabilizes 
through  collisions  with  the  environment.  We  will  indicate  this  by  “00”  as  superscript  {k°°{T)). 
For  the  formation  of  NO2  the  same  criterion  as  in  Ref.  [12]  is  assumed,  that  is,  the  oxygen 
atom  02  is  captured  once  it  approaches  to  within  less  than  3.78  uq  of  the  center  of  mass  of 
NOT 

(II)  is  for  NO2  formation  and  survival  which  carries  the  label  “NO2”  as  superscript 

For  this  case,  the  “oxygen  capture  criterion”  (see  item  (I))  is  combined  with  a  lifetime  crite¬ 
rion,  defined  as  tufe  =  Uraj  —  ttof  >  0.  Here,  ttraj  is  the  total  time  of  the  trajectory,  which  is 
either  (a)  ttraj  =  100  ps  for  the  case  that  NO2  does  not  dissociate  or  (b)  the  time  until  any 
of  the  interatomic  separations  is  larger  than  the  initial  separation  of  NO  1-02.  The  interac¬ 
tionless  time  of  flight  ttof  =  2y^ —  b‘^)/vmi  is  determined  from  the  initial  conditions. 

Case  (III)  is  the  oxygen  atom  exchange  reaction,  “ex”  with  associated  rate  coefficient  k^^(T). 
The  criterion  for  this  case  is  that  the  N-02  distance  is  small  and  the  N-01  or  01-02  distances 
are  large. 

Process  (IV)  refers  to  oxygen  molecule  formation  (0102-l-N)  and  the  corresponding  rate 
coefficient  is  {T).  This  case  is  encountered  if  the  01-02  distance  is  small  and  the  N-01 

or  N-02  distances  are  large. 


3.4  Monte  Carlo  calculation  of  the  thermal  rate  coefficients 

The  thermal  rate  coefficient  [23,  24,  25]  can  be  determined  from 


/c(re,Tt,r„) 


A  [m 

g{Te)\l  TTfl 


a{Ec]Tr^)Ece  ^^^‘^dEc, 


(5) 


where  /3t  =  k^Tt  and  ks  is  the  Boltzmann  constant.  Ft,  Tg,  T^v  are  the  translational  temper¬ 
ature  of  NO  and  O,  the  electronic  temperature  of  NO  and  the  rovibrational  temperature  of 
NO,  respectively,  g{T^)  is  the  electronic  degeneracy  factor, [11,  26]  g  is  the  reduced  mass  of 
NOl  and  02,  and  cj(Fc;Trv)  is  the  integral  cross  section  as  a  function  of  the  collision  energy, 
Ec  and  Trv.  In  our  studies  it  was  assumed  that  Tt  =  Te  =  T^.  In  other  words,  the  various 
degrees  of  freedom  in  the  reactants  are  in  thermal  equilibrium.  If  no  subscript  is  shown  for 
T  then  it  corresponds  to  the  common  temperature.  In  a  similar  manner,  the  reaction  cross 
section  a{Ec‘,Trv)  for  a  given  collision  energy  Ec  is 


a{Ec;Tr^) 


E 


,ax  +  l)e-^-^^^arj{Ec;  vj) 


(6) 
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where  a^jiBc)  is  the  state  dependent  cross  section  at  collision  energy  E^.  The  energy 

E^j  of  a  rovibrational  state  (f,j)  is  calculated  according  to  a  Morse  oscillator  model  for  the 
NO  molecule.  [22]  The  cross  section  as  an  integral  of  the  opacity  function  Pyj  (6;  Ec)  for  given 
Ec  and  rovibrational  state  (u,j)  is 


ayj{Ec)  =  J  Pyj{b-,Ec)2Trbdb  (7) 

0 

Assuming  a  thermal  distribution  of  the  reactants  0(^P)  and  NO(^n)  the  expression  for  the 
electronic  degeneracy  factor  g{Te)  is 


g{Te)  =  [5  +  3exp(— 227.8K/Te)  +  exp(— 326.6 K/Te)]  [1  +  exp(— 177.1  K/Te)].  (8) 

The  integral  in  Eq.  5  can  be  calculated  by  using  an  Importance  Sampling  Monte  Carlo 
scheme. [27]  For  this,  the  vibrational  and  rotational  quantum  numbers  v  and  j,  and  the  colli¬ 
sion  energy  (Ec)  are  sampled  according  to  Boltzmann  distribution. 

4  Results  and  discussion 

4.1  The  ground  state  PES 

The  smoothness  of  the  RKHS  interpolation  is  evident  in  the  contour  plots  of  Figure  2.  In 
panel  (a),  the  PES  for  the  “NOl-l-02”  channel  is  shown.  The  global  minimum  was  found  at 
i?i  =  3.21  ao  and  ai  =  150.8°  (rNO  =  2.25  uq  and  Oq-n-o  =  133.1°  (Table  1)  in  internal 
coordinates).  The  PES  for  the  “0102-l-N”  channel  is  reported  in  panel  (b).  Two  symmetric 
minima  are  observed  for  the  same  Jacobi  distance  which  is  measured  from  the  center  of 
mass  of  O2  to  the  N  atom. 

Some  features  of  our  ground  state  PES  are  summarized  in  Table  1  and  compared  to  re¬ 
sults  from  previous  calculations[13]  (using  icMRCIJ-Q/cc-pVQZ  with  CASSCF(13,10))  and 
to  experiments [28,  29,  30,  31,  32]  (high-resolution  microwave  spectroscopy[29]  and  laser  in¬ 
duced  fluorescence[31]).  Good  agreement  with  previous  computations  and  experiments  was 
found  from  the  methods  used  in  the  present  work. 

4.2  Thermal  rate  coefficients 

A  total  of  10  000  individual  trajectories  was  run  at  each  temperature,  including  T  =  100,  200, 
300,  400,  1000,  1500,  2000,  2700,  5000,  7000,  10000,  15000  and  20000  K.  In  the  following,  the 
thermal  rate  coefficients  obtained  from  the  trajectory  calculations  are  presented. 

In  Figure  3,  the  thermal  rate  coefficients  for  NO2  formation  in  the  high-pressure  limit, 
k°°{T),  are  shown  and  compared  with  published  computational  and  experimental  results. 
Even  though  the  major  interest  here  concerns  the  high-temperature  range,  these  data  provide 
a  validation  of  the  computational  procedure  put  forward  in  the  present  work.  We  report  the 
error  bar  as  twice  the  standard  deviation  (as)  in  order  to  include  95.4  %  of  the  cases  in  the 
confidence  interval.  In  the  following  it  is  important  to  emphasize  that  the  present  results  are 
from  simulations  using  the  l^A'  PES  only  (open  squares  with  solid  line).  The  present  results 
are  in  reasonable  agreement  with  results  from  previous  calculations  [11]  (solid  down-oriented 
triangles)  which  also  used  the  l^A'  PES  only,  although  with  a  different  parametrization. 
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Figure  2:  Contour  plots  for  the  NOl+02  (panel  a)  and  020H-N  (panel  b)  grids  in  Ri,  ai 
and  i?3,  as  Jacobi  coordinates,  respectively.  The  N-01  (ri)  and  02-01  (ra)  distances  were 
set  to  2.25  uq.  In  both  cases  the  contour  levels  are  separated  by  0.5  eV. 

Both  calculations  underestimate  the  experimentally  determined  thermal  rate  coefficients  at 
T  =  200,  300  and  400  K  (solid  circles  and  squares).  At  2700  K  the  computed  rate  coefficient 
compares  favourably  with  experimental  data.  [33,  34] 

No  correction  was  made  for  ZPE  and  other  quantum  effects,  such  as  tunneling  through 
the  centrifugal  barrier  which  can  be  still  significant  at  low  temperatures.  The  centrifugal 
barrier  originates  from  the  conservation  of  the  total  angular  momentum  and  consists  of  the 
rotational  energy  associated  with  the  orbital  angular  momentum  of  the  colliding  partner  with 
respect  to  the  center  of  mass  of  the  whole  system.  The  orbital  angular  momentum  is  directly 
proportional  to  the  impact  parameter  and  with  the  initial  relative  momentum  of  reactants. 
An  increase  of  the  thermal  rate  coefficient  for  lower  temperatures  has  been  found  to  be  related 
to  tunneling  through  the  entrance  barrier  in  the  reaction  of  OH  with  methanol  (T  <  200  K) 
and  0(^P)  with  alkenes  (T  r\j  20  K).[35,  36]  This  has  also  been  shown  in  theoretical  studies 
of  HCOH  isomerizations.[37]  On  the  other  hand,  results  when  using  the  two  lowest  l^A  and 
l^A  PESs  indicate  that  the  l^A  PES  might  play  an  important  role  along  with  quantum 
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Table  1:  Equilibrium  geometries  and  normal  mode  frequencies  from  the  present  work  and 
from  reference[13]  and  experiments.  The  variables  ri,e  =  r2,e,  Oo-N-o  and  De  are  the 
equilibrium  N-0  distance,  the  O-N-0  angle  and  the  well  depth  for  the  global  minimum  of 
NO2,  respectively. 


Equilibrium  NO2 

Present  work 

Ref.  [13] 

Experiments 

ri,e  =  r2,e  (ao) 

2.25 

2.26 

2.255  [29] 

0 

1 

1 

0 

0 

133.1 

134.3 

133.9  [29] 

De(eV) 

3.21 

3.11 

3.23  [30,  31,  32] 

Normal  frequencies  (cm 

Symmetric  stretch 

1322.7 

1333.6 

1319.794  [28] 

Asymmetric  stretch 

1623.3 

1658.0 

1616.852  [28] 

Bending  mode 

759.6 

755.0 

749.649  [28] 

effects  at  temperatures  below  400  K.  However,  for  temperatures  above  1000  K  the  results 
using  only  the  l^A  PES  reasonably  match  with  the  experimental  results,  e.g.  for  1500  and 
2700  K.  This  can  be  explained  by  noting  that  the  well  of  the  electronically  excited  C-state  of 
NO2  (C  l^A  )  is  comparatively  small  and  leads  to  redissociation  of  NO2  at  high  temperatures 
[11].  Given  this,  the  influence  of  tunneling  and  the  l^A"  PES  at  high  temperature  appear  to 
be  negligible. 

Eor  practical  work  in  reaction  networks  [39]  it  is  advantageous  to  use  a  parametrized  form 
for  k°°{T).  A  previously  considered  empirical  expression[ll]  has  been  fitted  to  the  results 
from  simulations  using  the  l^A'  PES, 

A:“(r)/(10-^^cmVs)  =  2.715  +  0.4490 InT/K  +  0.3385  {InT/Kf  -  0.7134  (InT/K)^  .  (9) 
This  yields  the  red  solid  line  in  Eigure  3. 

The  thermal  rate  coefficients  for  temperatures  between  T  =  5000  K  and  T  =  20000  K  are 
reported  in  Eigure  4  for  cases  I  to  IV  studied  in  this  work.  With  increasing  temperature  a  slight 
increase  in  the  thermal  rate  coefficients  is  found.  At  20000  K,  they  increase  to  =  5.5 T0“^^ 
cm^/s,  =  3.4- 10“^^  cm^/s,  =  1.9- 10“^^  cm^/s  and  =  0.65- 10“^^  cm^/s.  As 

could  be  expected,  results  with  and  without  ZPE  correction  agree  to  within  the  error  bars. 
Again,  for  practical  purposes  the  data  for  k°°{T)  were  fitted  to  a  linear  equation  as  (solid  line 
in  Eigure  4) 

A:“(r)  /(10-i^cmVs)  =  2.497  +  0.1529(r/K).  (10) 

Eor  a  more  detailed  description  of  the  O2  production  channel,  Eigure  5  (panel  (a))  re¬ 
ports  the  fraction  of  initial  vibrational  energy  of  NO  relative  to  the  total  energy  versus  the 
vibrational  quantum  number  of  NO  (unq)  for  8  values  of  the  collision  energy.  The  fraction 
of  collision  energy  is  complementary  to  the  vibrational  energy  to  give  the  total  energy  of  the 
system  {Ec  +  Evib  =  Etot)-  In  panel  (b),  the  contour  plot  of  the  channel  probability  is  shown. 
In  all  cases,  the  rotational  ground  state  (jno  =  0)  is  assumed.  Eor  this  calculation,  2500 
trajectories  were  run  for  each  combination  of  Ec  ={0.82,  1.09,  1.36,  2.18,  2.72,  4.08,  4.63, 
5.44}  eV  and  uno  ={0,  1,  2,  3,  4,  6,  9,  12,  15,  20}  values.  The  results  indicate  that  for 
=  0  the  O2— channel  opens  only  for  very  high  collision  energies  (>  4.63  eV).  However, 
for  an  initial  NO-state  with  uno  =  9  (corresponding  to  Ec  =  0.82  eV)  production  of  0.6  % 


Distribution  A:  Approved  for  pubiic  reiease;  distribution  is  uniimited. 


Figure  3:  Thermal  rate  coefficients  for  the  NO2  formation  in  the  high-pressure  limit,  k°°{T) 
versus  temperature.  In  the  legend,  [a]=Ref.[ll],  [b]=Ref.[12],  [c]=Ref.[10],  [d]  =Ref.  [33] , 
[e]=Ref.[34],  and  [f]=Ref.[38] 


of  O2  is  observed.  This  qualitatively  follows  from  Polanyi’s  rule  [23]  which  predicts  that  the 
reaction  should  go  to  the  late-barrier  “O2  +  N”  channel  when  the  NO  molecule  is  initially 
vibrationally  highly  excited  even  for  low  collision  energies.  In  order  to  give  a  clearer  picture 
of  this,  let  us  assume  we  start  the  reaction  coming  from  N-I-O2  (left  side  in  Figure  1).  In  this 
situation,  the  reactants  experience  an  early  barrier  of  about  1.3  eV  (Figure  1  and  3  b).  The 
exit  channel  in  this  case  is  barrierless  (NOl  -|-  02).  According  to  Polanyi’s  rule,  the  excess 
energy  released  after  crossing  the  early  barrier  goes  predominantly  into  the  internal  degrees 
of  freedom  in  the  products  (NOl)  rather  than  into  their  (relative)  translational  energy. 

For  the  set  of  vibrational  quantum  numbers  and  energies  studied,  the  maximum  channel 
probability  of  13  %  for  O2— formation  is  reached  at  uno  =  20.  In  general,  the  trend  indicates 
that  for  lower  vibrational  quantum  numbers,  more  energy  needs  to  be  in  the  translational 
degree  of  freedom.  Conversely,  for  higher  vibrational  excitation  of  the  NO,  less  collision  energy 
is  required.  When  both  the  collision  energy  and  the  vibrational  quantum  number  of  NO  are 
high,  the  reaction  probability  of  N-I-O2  formation  is  maximal. 

5  Conclusions 

The  0(^P)  -|-  NO(^n)  reaction  has  been  studied  over  a  range  of  temperatures  relevant  for  the 
hypersonic  flight  regime  (5000-20000  K).  For  this  purpose  the  thermal  cross  sections  and  rate 
coefficients  were  calculated  for  the  different  exit  channels.  The  results  indicate  that  oxygen 
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Figure  4:  High-pressure  limit  thermal  rate  coefficients  for  NO2  formation  (squares)  (fit  with 
black  solid  line),  NO2  formation  and  dissociation  (circles),  oxygen  exchange  channel  (up- 
oriented  triangles)  and  oxygen  production  channel  (diamonds). 


exchange  (N02-I-01)  occurs  at  all  temperatures.  However  O2— formation  is  only  found  for  T  > 
5000  K.  For  both  channels,  the  thermal  rate  coefficient  increases  as  the  temperature  increases 
(T  =  5000  to  20000  K).  A  more  detailed  scrutiny  of  the  “0102-l-N”  channel  indicates  that  the 
smaller  the  vibrational  quantum  number  of  NO,  the  higher  translational  energy  is  required  to 
open  it;  and  the  contrary,  the  higher  the  vibrational  quantum  number,  the  less  translational 
energy  is  needed.  Our  results  for  lower  temperatures  are  discussed  and  compared  to  previous 
experimental  and  computational  results.  In  this  respect,  this  work  also  presents  the  validation 
of  an  efficient  computational  protocol  for  the  investigation  of  reaction  dynamics  of  triatomic 
systems.  Two  aspects  make  the  approach  efficient  and  accurate.  Firstly,  the  ab  initio  energies 
are  handled  by  means  of  a  RKHS  representation  which  is  accurate,  computationally  convenient 
and  correctly  describes  the  asymptotic  regions  of  the  PESs.  Secondly,  thermal  quantities  are 
calculated  by  using  Importance  Sampling  which  is  efficient  to  evaluate  multi-dimensional 
integrals  over  the  phase-space  of  initial  states.  With  less  than  1000  classical  trajectories  a 
relative  error  below  10  %  for  the  rate  coefficients  is  obtained. 

6  Outlook 

The  present  work  is  the  basis  for  computing  state-selected  cross  sections  and  rate  constants 
for  the  reaction  considered.  More  broadly,  the  computational  framework  discussed  here  can 
be  generalized  to  other  relevant  reactions  for  the  hypersonic  flight  regime.  Working  along 
the  lines  put  forward  in  the  present  report,  starting  from  electronic  structure  calculations, 
building  the  fully-dimensional,  reactive  PES  is  almost  automatic  as  is  generation  of  the  initial 
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Figure  5:  Fraction  of  initial  vibrational  energy  of  NO  with  respect  to  the  total  energy  versus 
the  vibrational  quantum  number  for  8  values  of  the  collision  energy  (panel  (a)).  The  contour 
plot  of  the  channel  probability  as  percentage  is  shown  in  panel  (b).  Only  data  for  “O2+N” 
channel  are  shown. 


conditions  and  running  the  necessary  atomistic  simulations.  This  provides  the  raw  data  for  the 
computation  of  the  rate  coefficients  at  temperatures  for  which  concrete  experiments  is  difficult 
and  which  is  the  major  information  sought  for  incorporation  into  the  chemical  modeling. 
Pursuing  this  for  additional  relevant  reactions  and  including  the  effect  of  electronically  excited 
states  will  be  of  particular  interest  for  the  continuation  of  this  project. 
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